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Within the modified model of proton ordering of one-dimensional ferroelectric having hydrogen bonds of 
PbHP04-type, their thermodynamic and dynamic characteristics are studied and calculated taking into account 
the linear (by crystal deformations c/ (i = 1,3) and £4) contributions into the energy of a proton system but 
without taking into account the tunneling in the two-particle cluster approximation. There has been obtained a 
good quantitative description of the temperature dependence of polarization, static dielectric permittivity, heat 
capacity and frequency dependence of dynamic dielectric permittivity at different temperatures for PbHP04 
and PbHD04 crystals. 
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1. Introduction 

Ferroelectric properties in PbHP 04 (LHP) and PbDP 04 (LDP) crystals were disclosed in work jll]. At 
temperature - 310 K in LHP and T®' = 452 K in LDP a phase transition of the second order takes 
phase. These are monoclinic crystals of P2/c space group in paraelectric phase jUil . Ferroelectric phase in 
LHP is characterized by a Pc symmetry having a spontaneous polarization in the direction that forms an 
angle of 10° with the crystal a-axis. The elementary cell of LHP contains two molecules. The parameters of 
unit cell of LHP are as follows a = 4.688 A, b = 6.649 A, c = 5.781 A, ;S = 97.11° and of LDP — a = 4.6855 A, 
b = 6.6911 A, c = 5.7867 A, ;S = 97.10°. 

A characteristic feature of crystal structure of ferroelectric of LHP is the presence of hydrogen bonds 
that link the PO 4 tetrahedrons into infinite chains that stretch out along the c-axis. According to the 
two possible equilibrium positions, the protons (deuterons) at these bonds in a paraelectric phase are 
distributed statistically uniformly, whUe in ferroelectric phase there appears a spontaneous asymmetry 
of population. The data presented in works dSl testify to the fact that the proton on the hydrogen O- 
H... O bonds move in two-minimum potentials. A noticeable change of the phase transition temperature 
at deuteration of LHP |ll-@], as well as of dielectric (301 and thermal @] properties testifies to an impor¬ 
tant role of a cooperative behavior of protons in the appearance of ferroelectricity in this type of crystals. 

Taking into consideration that the direction of dipole moment of V 4 {Bu) internal mode is close to 
the direction of spontaneous polarization Ps i^[l 3 l. as well as since this mode shows a considerable 
temperature softening in a ferroelectric phase, which is similar to the temperature dependence Ps{T), we 
may assume that this internal mode V 4 plays a crucial role in the mechanism of the phase transition in 
LHP crystal. The inter-phonon interaction between the low-frequency Ws and internal V 4 modes that have 
got the same symmetry (Bu), close to Tc, may result in static deformation of PO 4 groups, which causes 
a static dipole moment of V 4 mode. It is natural to assume that the value of Ps is mainly determined by 
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a static dipole moment of PO4 jroups. The angle between these moments and the vector of spontaneous 
polarization is 2°. The work (IZD estimates the contribution of the polar deformation of PO4 groups into 
LHP to be ~ 80% of the experimentally observed value of LHP. The phase transition into LHP takes place 
due to the proton ordering, phonon anharmonicity and proton-phonon interaction. 

Microscopic models of phase transition of LHP were discussed and studied in works The 

model discussed in iish corresponds to LDP. Works present the LHP model that takes the 

tunneling of protons on hydrogen bonds into consideration. Work 12111 considers a simple two-sublattice 
model of partially deuterated crystal. In should be noted that works Im|2i 1] take proton-lattice interaction 
into consideration as well. Moreover, work (3 considers the anharmonicity of the lattice oscillations. 
Unfortunately, in [3. there was used an approximation of mean field while studying the static 
and dynamic dielectric permittivity, which is insufficient for LHP. In work fii. a simple model LHP that 
takes tunneling into consideration is solved within the approximation of two-particle cluster. Using the 
Green’s function method, the work fl^ presents temperature dependencies of relaxation times of LHP 
and LDP, while work El presents dynamic dielectric permittivities. However, all these works do not put 
forward the task of describing the corresponding experimental data. 

The model of a deformed crystal was beneficially used for a description of dielectric, piezoelectric, 
elastic, thermal and dynamic characteristics of quasi-one-dimensional ferroelectric with hydrogen bonds 
of CSHPO 4 type in (z^- 

This work presents a modified model of proton ordering of hydrogen bonded one-dimensional ferro¬ 
electric of LHP-type that takes into consideration the linear [by deformations £,■ {/ = 1, 3), £ 4 ] contribution 
into the energy of the proton system. Within the approximation of two-particle cluster, their dielectric, 
piezoelectric, elastic, thermal and dynamic characteristics are calculated. 


2. Model Hamiltonian of PbHP 04 crystal 


Let us consider a proton subsystem LHP that moves on the 0-H... 0 bonds forming zigzag-like chains 
along the crystal c-axis. 

The unit cell of LHP is formed by the chain that contains two neighboring PO 4 tetrahedrons together 
with two short hydrogen bonds that relate to one of them (tetrahedron of ‘A’ type) (see figure [1). The 
hydrogen bonds that adjoin the second tetrahedron (of ‘B’ type) belong to two closest structural elements 
of A’ t)rpe that surround it. 

The Hamiltonian of the proton-ion system of LHP, neglecting the tunneling effects of the protons on 
the 0-H... O bonds, is as follows izzilz^l : 

H = lVr’(/iseed-2mi^-^-^(dR^R ,-i-5R^+r,R Y, 

+ Y(>)dk]b'libici + ^ EE E T//(fc)e‘“'' [bu + + 

kl V iV qf kln=l V Z Z ; 

“EEfM/(fc)cos5°£’4 + ^ii{k)cos2°Es] (hfc/ -1- (2.1) 
q kl^ ' VN 


where v is the volume of unit cell, N is the total number of unit cells, cr^f is the z-component of pseudo¬ 
spin operator that corresponds to the proton located in the ^ 7 -th cell on the /-th bond (/ = 1,2). Eigenval¬ 
ues of the operator aqf = +1 correspond to two possible positions of the proton on the hydrogen bond. 
Parameters /ip and /i/ (fc) are dipole moments that correspond to the proton on the hydrogen bond and to 
the dipole-active phonon mode, respectively. Uiseed is the seed energy that appears in the form of crystal 
deformations £/, £4 and electrical field along the crystallographic axis a* which is perpendicular to 
the plane (b,c), respectively, and contains elastic, piezoelectric and dielectric parts: 


1 3 

f^lseed - T Y 
^ i,i=l 


/=! 


.£0 


M + -ci4°£4 - Y - e?4^l£4 “ 

^ i = l 


( 2 . 2 ) 


where cf9, cf°, cf®, e°., are seed elastic stress, coefficients of piezoelectric strain and dielectric 

susceptibility of a mechanically clamped crystal. 
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Figure 1. Unit cell of LHP crystal lisll . 


The second term in I I2.H is the Hamiltonian 
of short-range interactions between protons. The 
first Kronecker’s symbol corresponds to the pro¬ 
ton interaction in the chain close to the tetrahe¬ 
dron of ‘A’-type, the second Kronecker’s symbol 
corresponds to proton interactions near the ‘B’- 
type tetrahedron, r is radius vector of the relative 
position of the proton bond in the unit cell. The 
quantity mi that describes short-range proton in¬ 
teractions within the chain may be expanded into 
a series with respect to deformations e;, £4 re¬ 
stricted to the linear summands jj^ : 

3 

Wi - W+J^Si£i + S 4 £ 4 . (2.3) 

1=1 


The third term describes the effective long-range 
dipole-dipole interaction between protons within 
the chain running along c-axis, the fourth summand represents the lattice vibration energy b^i — 
Bose operators, /-index of phonon branch), the fifth summand corresponds to the proton-phonon inter¬ 
actions. The last term describes the interaction of the lattice with the external electrical held. 

Long-range interactions of protons and interactions of protons with the lattice vibrations are consid¬ 
ered within the the mean field approximation. Thus, the Hamiltonian (2.1) looks as follows: 


H=Hp + Hi, 


(2.4) 


where 


Hu = NvUisee<i- 2 wiY,-^-^[SRqR^HSR^+r.R^,] 


qq‘ 


1 (CT'”h (CT%) (a'”> ) , 

+ 0 L -^- L ^-^ 

^ qq'.ff ^ qq'Jf ^ 

^ \ if 


y/N 


^ ^ T//(fc) (h fc;-H J.,) - 


(2.3) 


qf kl 


kl vNqfkl 

-'Z.'Z.\l^lik)cos5°El -i-p;(fc)cos2°£’J {bki +b^^i). (2.6) 

q kH ' yN 

Using the Heisenberg equation of motion for the mean values of Bose-operators 

^,(bti)^-i{lbt,,H]), (2.7) 


we find that 


(bki + b^kl> 


_L y ^-ika„,'^q'f' 

q'f d)l(k) 2 

-I- y' - 1 ;- [/i;(fc)cos 5°£'i* + iii(k)cos 2°Eg] —— . 

^ J ^/N 


( 2 . 8 ) 
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Taking into consideration the expression <2.5> . the Hamiltonian of system Hp takes the following form: 

D . + (5lJ _i_^ D .1 

n' 2 2 


Hp = NvUiseed-2wiY,-^-^{SRqR^:+^Rq+r.R^,) 

qq’ 


+ ^ E - y -E -T 

qq'Jf ^ qq'.ff ^ 

-^[^cos 5° El +HCOS 2°Es]i-^ + -^Y 

cj '• Z Z •' 


The following notations are used herein: 


Jffiqq'] 

E 


1 ^ 2Tif{k)Tlf{-k) 


Jfriqq'] + —Y. 

q’ kl oiiik) y/N 


Considering the symmetry of unary function of proton distribution 


= (CTql) = ((Jq2} 


(2.9) 

( 2 . 10 ) 


and expanding the constant of long-range proton-proton interactions into a series with respect to defor¬ 
mations £;, £ 4 , restricted to the hnear summands: 


3 


}ff'^}+Y,Vu£i + Vi£4, 

i=l 


( 2 . 11 ) 


where Jfp = Zr^-r^, Jffiqq') are Fourier transforms of the constants of long-range interaction, we get 
an output Hamiltonian H in the following form: 


H — A(nf7iseed 


O'ql CTq2 


+ Sr^+t.r^,) - 2Nvi (j7 '1’^)2 - -t 


qq 


aq2' 


( 2 . 12 ) 


q ' ^ ^ q 

where the following notations are used: 


Vl = V-t^ 1 /^,£/-|- 1 // 4 £ 4 , 


/ij = ficos 5°, /is = /icos 2°. 


The approximation of two-particle cluster is used in order to calculate the physical characteristics of 
PbHP 04 -type compound. In this approximation, the LHP thermodynamic potential is as follows: 

Gi = NUueed + 2Nvi -ksTj^l )nSpe~^^4A + InSpe”^"^® - - InSpe”^"™ 

q I 2 

-Nv’^aiEi - Nva^Ei, (2.13) 

i 

where H^q\f^gy Hq^ are two-particle and one-particle Hamiltonians assigned by the following equations: 




- CT™ 


^ (Jq2 \ 

P y 2 2 I’ 


(2.14) 

(2.13) 
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The following notations are used herein: 

XE^P[-A‘' + 2 viT]^^'>^ + nlEl+HsEs), xe^-PA'^ + xe, ; 6 =l/A:Br, (2.16) 


where A“ is the effective field formed dy the neighboring links beyond the cluster borders. Within the 
cluster approximation, the field is determined according to the condition of self-consistency, i.e., the 
mean value of pseudospin should not depend on a particular Gibbs distribution (two-particle or 

one-particle Hamiltonian) according to which it is estimated: 


Spa®e-Ki Spr7»>,e-K 


qf 


Spe 




qA 


Spe 




(2.17) 


Then, based on <2.17t and considering (2.141 and (2.151 . we get an equation for the mean value of 
pseudospin in the following form: 




sinhx^ sinhx^ 


ai + coshxg 


D 


(2.18) 


where 


1 . l + , ( 11 , 

«' 2'-w' ’ + +— 


ai = exp 


-P 


3 

w+ '^Siei + 64e4 


i=l 


3. Static dielectric, piezoelectric, elastic and thermal characteristics 
ofPbHP 04 


Having calculated the eigenvalues of two-particle and one-particle Hamiltonians, let us present the 
thermodynamic potential (2.131 per unit cell in the following form: 

Gi 

gi = — = n(/iseed-2fcBrin2-(m-i-5iei-i-52e2 + 53e3 + 55£5) 

+ {v + y/iEi + y/2E2 + V3£3 + y/ 4 £i]iri^^^^f 

- ^bTIh [l - (p^^^^)^] - 2 A;Brin(ai -i- coshx^) - v((JiEi + (T 2£2 + ce^es] - va 4 E 4 . 


Using the equilibrium equation 


vidEiJEi v[dEi 


we get an equation for deformation e/, £4 and polarization Pi: 


n — , £0 , £0 , £0 c _ p 

0 - £i + C -2 £2 + ^3 + ^14 ^4 


n — , £0 , £0 . £0 /^O p 

0 — Ci4 £i + C24 £2 + C34 £3 + C44 £4 — ei4ti — 




1 [ 5gl ] 
v[dE;j^^ 

-PI, 

Si coshx- ai 

V 

V ai + coshx 
S 4 coshx - a\ 

V 

V a\+ coshx 


(3.1) 

(3.2) 

(3.3) 

(3.4) 


Based on the relation (3.11 and (3.41 . we get the following thermodynamic characteristics of LHP crys¬ 
tals, i.e., isothermal static susceptibility of a mechanically clamped crystal: 


xlf- 


dP* 
OF* , 

"-^1 'si 


Pxe 

■ All + ■ 


2 v D-xe<P^ 


(3.5) 
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where the following notations are used: 

xe - coshx£ - 77^^'^sinhx£, 
isothermal coefficients of piezoelectric stress: 


<P 


n - 




+ Pvi; 


^ piy/jx + djai) pT _^o + 


[dsilEi ^ D-xip^ 

isothermal elastic constants at a constant field: 


■ 


de4 j 


El 


^(1). 


V D- xqf! 


dCFj 


cl 


EO s s «lC0shx P T] 


( 1)2 


^cfr-^SiSj 


i I El 


D2 


V D-xqf! 


SiSj —^ + y/iy/jX + (Siy/j +6 jy/ijai 


E _ EO )S<5i54COShx Pt]' 


( 1)2 


vD^ 


EO _ EO PS^coshx 


v[D-x(p^) 

p^m2 


a\q)^ 


6i8i + ^fiipiX + [Siy/4 + y/i64]ai 


vD^ v[D-x(p’^) 


, ^,,2 


D 


+ ^ij^x + 254i/i4ai 


Other dielectric, piezoelectric and elastic characteristics of LHP may be calculated using the above 
mentioned results. In particular, we may get the matrix of isothermal compliance at a static field sf., 
which is inverse to the matrix of elastic constant cf .; isothermal coefficients of the piezoelectric strain: 


4 = Z44- (U=1,2,3,4), 


(3.6) 


isothermal dielectric susceptibility of free crystals: 

i 4 = ifr+E 44 - 


(3.7) 


In the LHP crystals, there takes phase a phase transition of the second order from paraelectric phase 
into ferroelectric phase at temperature that satisfies the equation: 


f mi ) vi 


(3.8) 


The molar entropy of a crystal conditioned by a proton subsystem is obtained through a direct differen¬ 
tiation of thermodynamic potential: 


S 


■R [ dgi) 
2 


^ I In 2-1 In [ 1 - (qd) )2 ] + 2 In D - 2j0Vi (qd) )2 + I ^ 


( 3 . 9 ) 


where f? is a universal gas constant. 

The molar heat capacity of a hydrogen subsystem of LHP at a constant strain is calculated by a direct 
differentiation of entropy d3.9L 


ACp 


T 

T 





-1 T 






dS ) I dEi ) 

de4J^idr )' 


(3.10) 


4. Relaxational dynamics of PbHP 04 crystal 

This section describes the dynamic phenomena in LHP at the application of electrical field to a 
crystal. While calculating the dynamic characteristic of this kind of ferroelectrics we use a kinetic equa¬ 
tion based on the method of non-equilibrium statistical Zubarev operator liill. 
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The kinetic equation for the mean values of pseudospin operator is as follows: 


, <Pm> - 


Qqftia ^Pm) + tanh— [pm) 


k;, 


where 


C50 

J<C“ = J dte~^^cosD.‘^tRe{u[t)u'^}q, o; = 0 ,+l, 


(4.1) 

(4.2) 

(4.3) 


while {u“{t)u“')q are correlation function of a thermostat, is a Fourier component of the op 


erator a“At), Q“ are eigenfrequencies of the Hamiltonian of pseudo-spin model 1 12.14) . o"'', = aqf, 


qf 

+ r ■ V 

O'-„ = a f + ia 

qf qf qf 


Taking into account the time evolution law of pseudo-spin operators and using the frequency 

presentation of these operators, one may calculate the expressions for Q^f^a^Pqf) iii another form. Fi¬ 
nally, it allows us to rewrite the kinetic equation )4.1) as follows: 


d 

dt 

d 

dt 


-t + bi, 

- 77 '^)^ = bzipP^^ -I- + bz, 

as well as the equation for a unary function in one-frequency approximation: 


— 77 ™^ = - 2 ^ 077 '!’^ -H 2 )<:otanh— 
dt 2 


(4.4) 


(4.3) 


The following notations are used herein: 

hii = -(^Ti + K-i) + K1Z1- K-iZ-i, 

bzi = - 2 (^:i - ^:_i) + 2(^1 Zi + r_iZ_i), 


hi2 = -^:i + ^:_i, 
(722 =-2(^1+ r_i), 


bi^KiZi + K-iZ-i, 
bz = 2{KiZi-K-iZ-i), 


Kp,- + K^- J dte ^^cos(0^t)Re{(M {t)u^) q +{u^ [t)u )q\, Zp - tanh ^^ ^ . 

0 

In case Kq - K-i - Ki- 112a, the system equations obtained in this section are in agreement with the 
equations obtained within the framework of a stochastic Glauber model (s^. Glauber equations describe 
a physical situation at which the Fourier images of the thermostat correlators are independent of the 
frequency (^1^ . 

Thus, from equations )4.4) and ( 14.5) . we find that 


where the following notation is used: 


1 

P= - 
2 



= - (1 - + L, 

dt 

a— 77 *^'^ = - 77 ^^^^ -I- tanh—, 
df 2 


tanh 

'pm xe^ 

-tanh 

/ Pwi Xe\ 


[2 2 ; 


[2 2 jj 

tanh 

' pWl Xe' 

-i-tanh 

pWi Xe\ 


,2 2 ) 


2 2 J. 


(4.6) 


(4.7) 
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Solving equations <4.61 in the case of small deviations from the equilibrium state, one may obtain the 
complex dielectric permittivity of the hydrogen subsystem of LHP: 




1 + (27rvT)2 ’ 




4nx*7ivT 
1 + (271VT)'^ 


(4.8) 


with the set of following notations 


X =■ 


fifP f2r[l-P™]-F 


2 v 


Yr 


-PiV + Vab) 


-1 


_i_ 1 Yr |2r[l-P™]-F 




a 2r-Y 


Yr 


-P[V+Vah)\, 




p(0) 

p(l) 


a\ 

1 + flj + 2flicoshjc ’ 
4ai(l- fll^)sinhjc 
[l + af + 2flicoshjc]2 ’ 


p(0) 


2flisinhx 
1 + + 2flicoshx ’ 

4fli [2fli + (1 + flj)coshjc] 
[1 + af + 2flicoshx]2 


5 . Comparison of numerical calculations with experimental data. 
Discussion of the results obtained 

Prior to the discussion of the developed theory, it should be noted that this theory, strictly speaking, 
holds for deuterated ferroelectric LDP. Thermodynamic and dynamic characteristics of hydrogen-bonded 
ferroelectrics taking tunneling Q into account, are essentially defined by an effective parameter of tunnel¬ 
ing O, which is renormalized by short-range interactions 13111 . Here, Q « n, i.e., an essential suppression 
of tunneling by short-range interactions takes place. Then, let us assume that the theory proposed by us 
holds for LHP crystals as well taking into consideration, in particular, the relaxational type of dispersion 
in LHP. 

Unfortunately, the elastic constants of the LHP crystal have not been experimentally determined so 
far. That is why it is impossible to specify the seed elastic constants and hence to calculate, based on 
the proposed theory, the piezoelectric coefficients, susceptibility of a mechanically free crystal, elastic 
constants in ferroelectric phase. In order to perform the numerical calculations of temperature and fre¬ 
quency dependencies of respective physical characteristics of LHP, the values of the following parameters 
should be specified: 

• parameter of a two-particle cluster w; 

• parameter of long-range interaction v; 

• effective dipole moment ps; 

• the seed dielectric susceptibilities Xu', 

• parameter a that defines the time scale of relaxation processes. 

In order to determine the above mentioned parameters, let us use temperature dependencies of ex¬ 
perimental physical characteristics, namely Tc, PsiT) d, eii(0,r) d . The value of the effective dipole 
moment is determined through the agreement of the theory with the experiment for polarization 
of saturation. In paraelectric phase, we determine by agreeing the theory with the experiment for 
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Table 1. A set of parameters of the theory for LHP and LDP crystals. 



Tc 

(K) 

w/ks 

(K) 

viks 

(K) 

Ms, 10-^“ 
(esu ■ cm) 

Mip. 10-^« 

(esu ■ cm) 

yCU 

All 

Ps 

(s) 

Rs 

(s/K) 

(s) 

Rp 

(s/K) 

PbHP 04 

310 

850 

19.98 

1.00 

1.46 

0.716 

0.339 

-0.009 

0.745 

0.002 

PbDP 04 

452 

1450 

2.19 

1.05 

2.00 

0.716 








Figure 2. Temperature dependencies of spontaneous polarization of LHP — o Ql.oQl and LDP —o Q], 



Figure 3. Temperature dependencies of dielectric permittivities of LHP — o Q], □ (^, a and LDP — 

A (li. □ (32I1. 



Parameter a is defined from the condition that 
theoretically calculated curves of should 

agree with the experimentally obtained curves. In 
is assumed that parameter a slightly changes with 
temperature: 


0 := [P + J?(Ar)]-10 


-14 


AT^T-Tr 


The unit ceU volume of LHP is taken to be equal 
to n = 0.1788 ■ 10“2i cm^, LDP — n = 0.1800 ■ 
10“2i cm^. 

The set of optimal parameters obtained this 
way is presented in table 1. 

Let us now discuss the obtained results. Fig- 
ure|2]shows temperature dependencies of sponta¬ 
neous polarization of LHP and LDP crystals together with the experimental data. It is seen in the figure 
that the data in papers Q] and 0] disagree between themselves. There is a good description of tempera¬ 
ture dependencies of spontaneous polarization obtained in paper |^. Polarization of saturation increases 
at the growth of the degree of deuteration x. 


Figure 4. Temperature dependence of heat capacity 
for LHP crystal: *, @|. 
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Figure 5. Frequency dependencies of real part of dielectric permittivity of LHP at AT: a) -10 K; b) 
-5 K; c) 5 K; d) 7.16 K; e) 10 K; f) 20 K and experimental data: □ — (3, o — (Hi, o — H, v — S; > — 

m. A — isl. 


Figure [3] presents temperature dependencies of static dielectric permittivities of mechanically 
clamped £^^(0, T) crystals of LHP and LDP calculated based on the microscopic theory, as well as the 
results of experimental studies 

As it is seen in figure the results of theoretical calculations of cJ'j (0, T) on the whole show good 



Figure 6. Frequency dependencies of imaginary e'j'j part of dielectric permittivity of LHP at A T: a) -10 K; 
b) -5 K; c) 5 K; d) 7.16 K; e) 10 K; f) 20 K and experimental data: □ — (sl, o — 0; ° — (1; V — El; > — 

El. A — E3i- 
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Figure 7. Temperature dependencies of real T) and imaginary ej\(v, T] parts of dielectric permit¬ 
tivity of LHP crystal at different frequencies v (GHz): 1 — 1, □ 10 — 2,D (^; 94 — 3; <> |3; 179 — 4, 

0 (Hi; 240 — 5, o 0; 480 — 6, o 0. 

quantitative agreement with experimental data obtained by the authors 

The temperature dependence of heat capacity for LHP crystal, together with experimental data of 
paper 0 are presented in figure|4] 

The dashed line shows the effective lattice contribution Co into the heat capacity which is estimated 
by us as the average of the difference CpCT) - ACpCT). Based on the proposed model and using the theory 
parameters (table [1), a quantitatively good description of the data of paper 0 has been achieved. The 
amount of the calculated heat capacity jump well correlates with the experiment. 

The measured temperature dependencies of a real and imaginary parts of dielectric permittivity at 
different frequencies for LHP are presented in papers ia l^l34l - [38ll . 

Figures!^ and [^present the results of calculations of frequency dependencies e[^(v] and r'/jlv), re¬ 
spectively, for LHP as well as experimental data. 

In these figures it is seen that there is a significant data scattering obtained in the experimental papers 
The best agreement for our calculated data of and and the experiment was 
reached with the results published in papers 

Figure [7] presents temperature dependencies of real e[^{v,T) and imaginary e' 1 ^(v,T] parts of dy¬ 
namic dielectric permittivity at different frequencies of LHP crystal as well as the experimental data. 

As seen in this figure, there is quite acceptable agreement between our theoretical results and the 
majority of experimental data. A significant disagreement is observed only with the experimental data of 
Briskot et al. (3^. The reason is a rather poor experimental accuracy (order of 40-30 %) of their experi¬ 
mental setup. 


6. Conclusions 

The present paper, based on a modified model of proton ordering that does not take into consideration 
the proton tunneling on hydrogen bonds in the approximation of a two-particle cluster, describes the 
theory of thermodynamic and dielectric, piezoelectric, elastic and dynamic properties of one-dimentional 
ferroelectrics of PbHP 04 -type. Optimal sets of model parameters have been found that make it possible 
to describe the available corresponding experimental data for LHP and LDP crystals. 
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TepMOAi/iHaMiHHi Ta Ai^HaMiHHi Aic/ieKTpi/iMHi B,nacTi/iBOCTi 
OAHOBI/IMipHI/IX CerHeTOeyieKTpi/IKiB 3 BOAHeBI/IMI/l 3B'B3KaMI/l 
Ti/iny PbHP 04 

I.P. BaneP, P.P. /leei/mbKi/iiP, 51.1/1. IHypP, O.B. BiReHbK^^ 

^ iHCTWTyr cJiSMKM KOHfleHCOBaHMX cncreM HAH yKpaiHi/i, syn. I CBeHqiqbKoro, 1,79011 /IbBiB, VKpaiHa 
^ HaqiOHanbHuit yHiBepcnrer "/IbBiBCbKa noairexHiKa", Bya. C. BaHflepw, 12, 79013 /IbBiB, VKpaiHa 

y paMKax MOflncJiKOBaHOi MOfleai nporoHHOro BnopaflKyBaHHa 3 BpaxyBaHHaM BsaeMOflii nporoHiB 3 HopMaab- 
HHMH KoaHBaHHBMn rpaTKia OflHOBWMipHnx cerHeroeaeKTpwKiB 3 BOflHeBHMH 3B'a3KaMM rnny PbHP 04 3 Bpaxy- 
BaHHnM aiHiitHHX 3a flecJopMaqiaMH Kpwcraay e; i £4 BHecxiB b eHepriK) nporoHHOi cncreMH, aae 6e3 BpaxyBa- 
HHa tyHeaiOBaHHn b Ha6an>KeHHi flBOHacrnHKOBoro Kaacrepa poapaxoBaHO i flocaiflxeHO ix repMOflHHaMiHHi 
i flHHaMiHHi xapaKTepncTHKia. OrpHMaHO flo6pnCi xiabKiCHWH onnc reMnepatypHOi aaaexHOcri noaapnaaqii, 
cranaHHOi flieaeKrpwHHOi npoHMKHOcri KpncraaiB PbHP 04 ra PbHD 04 , renaoeMHOcri i nacTOTHOi saaexHOcri 
flHHaMiHHOi flieaeKTpHHHOi npoHHKHOcri npn piSHiax reMnepatypax Kpncraay PbHPOa. 

K/noHOBi ctioBa: cerneToeneKTpi/iK, flieneKTpi/mna npoHiiKHicTb, n'eaoe/ieKrpi^i^Hi/iCi KoecpipieHT, PbHPO^ 
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